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Abstract 


This study aimed to investigate the occurrence of coronaviruses (CoVs) in captive birds placed inside a zoo- 
logical park in Brazil. The role of captive birds in the epidemiology of CoVs in the tropics is poorly understood. 
A total of 25 (1= 25) different species were tested for viral RNA using individual fecal samples collected from 
healthy birds. Reverse transcription—polymerase chain reaction targeting the 3’ untranslated region was used to 
detect CoV RNA, and positive samples were submitted for sequence analysis. The phylogenetic search revealed 
nine mutations in the black shouldered peafowl (Pavus cristatus) CoV sequence, which clustered separately from 
samples previously described in England. This is the first report on the detection of the CoV genome in captive 


birds in Brazil. 
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Introduction 


HE OUTBREAKS OF HIGHLY PATHOGENIC avian influenza 

viruses around the globe in the past few years have been 
devastating for the countries concerned and have created 
fear in the general population. Similarly, in 2002 the world 
witnessed the emergence of severe acute respiratory syn- 
drome in Guangdong Province, China, caused by a previ- 
ously unknown coronavirus (CoV) (Liu 2003). It is likely that 
additional, yet to be identified CoVs are circulating in do- 
mestic birds, especially in tropical countries. The three ex- 
isting groups of CoVs were initially distinguished based on 
the lack of an antigenic relationship between them (Masters 
2006). Group 3 CoVs, which are genetically similar to in- 
fectious bronchitis virus (IBV), have been detected in healthy 
galliform and nongalliform birds (Cavanagh et al. 2002). 
Recently, populations of wild birds in northern England 
were found to have IBV-like RNA in their feces (Hughes 
et al. 2009). Currently, little information is available about 
the role of South American bird populations as reservoirs or 
long-distance vectors of CoV-like viruses. 

Viral diseases can cause substantial mortality in both wild 
and captive bird populations. Brazil is a country with many 
migrating wild birds, which is likely to facilitate disease 
transmission from migratory to nonmigratory flocks in zoo- 
logical parks or on farms. As viral detection and identification 
techniques become increasingly sophisticated, and as the 


study of wildlife diseases increases, it is likely that new 
viruses will be discovered and new diseases will be encoun- 
tered (Hughes et al. 2009). Therefore, studies of the viral dis- 
eases of wild and captive birds are necessary. The latest 
techniques of molecular epidemiology can help scientists 
understand how viruses move from species to species and 
from region to region (Hughes et al. 2009, Jackwood et al. 
2010). 

The present study demonstrates the molecular detection 
and characterization of viruses from captive exotic bird pop- 
ulations sampled in Sao Paulo State Zoological Park in Brazil, 
which appear to be related to group 3 CoVs. 


Materials and Methods 
Study design 


Serial cross-sectional surveys of captive wild bird pop- 
ulations, maintained in a zoological park in the city of 
Bauru (N22°10’ 549°00" W), were undertaken in July 2009. 
Feces were collected from 25 different species of healthy 
birds (Table 1), placed in virus transport media (Eagle’s 
minimum essential medium supplemented with 10% fetal 
calf serum, penicillin, and streptomycin suspension, plus 
amphotericin B), and frozen at —80°C for individually 
analysis later. This work was authorized by the Animal 
Ethics Committee of the SAo Paulo State University (pro- 
tocol number 2008/09878-09). 


"Laboratory of Animal Virology, UNESP, School of Veterinary Medicine, University of SAo Paulo State, Aracatuba, Brazil. 
“Coronavirus Research Group, School of Veterinary Medicine, University of SAo Paulo, Sao Paulo, Brazil. 
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TABLE 1. CAPTIVE BIRD SPECIES FROM WHOM FECAL SAMPLES WERE COLLECTED AND SCREENED FOR CORONAVIRUS RNA 
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Common name Latin name Age Positive/total samples 
Chestnut-bellied heron Agamia agami Ad.* 1/20 
Egyptian goose Alopochen aegyptiacus Ad. 0/20 
Amazon parrot Amazona estiva Juv.” 0/5 
Blue-front parrot Amazona vinacea Ad. 0/5 
Bar-headed goose Anser indicus Juv 0/5 
Emperor penguin Aptenodytes patagonicus Ad. 1/2 
Cattle egret Bubulcus ibis Juv 0/0 
Red-legged seriema Cariama cristata Ad. 0/4 
Kelp gull Larus atlanticus Ad. 0/2 
Silver pheasant Lophura nycthemera Ad. 0/1 
Tawny-bellied screech-owl Megascops watsonit Ad. OF 12 
Red crested cardinal Paroaria coronata Juv 0/2 
Black shouldered peafowl Pavus cristatus Ad. 2/5 
Indo-Chinese Green peafowl Pavus muticus Ad. 0/14 
Ring necked pheasant Pavus torquatus Ad. 0/1 
Dusky legged guan Penelope obscura Ad. 0/1 
Rusted margined guan Penelope supercilliaris Ad. 1/2 
Pheasant Phasianus colchicus Ad. 0/12 
Greater flamingo Phoenicopterus rubber Ad. 0/9 
Black-fronted pippin-guan Pipile jacutinga Juv 0/8 
Caraca cheryway Plyborus plancus Ad. 0/23 
American purple gallinule Porhyrula matinica Ad. 23 
Saffron toucanet Pteroglossus bailloni Ad. 0/16 
Lettered aracari Pteropglossus inscriptus Juv. O12 
Toco toucan Ramphastos toco Ad 0/10 
Channel-billed toucan Ramphastos vitelinus Juv. 0/10 
Greater rhea Rhea americana Ad. 0/10 
Roadside hawks Rupornis magnirostris Ad. 1/10 


“Ad., adult (>2 years). 
Juv., juvenile (<2 years). 


Viral RNA extraction 


Two grams of feces was homogenized in two volumes of 
Eagle’s essential medium (Invitrogen™; pH 7.4) and clarified 
by centrifugation at 3000 g for 20 min. The supernatant was 
filtered once through a 0.75-um paper filter (Millipore™) and 
then twice through a 0.25-uwm syringe filter (Corning™) as 
described previously (Teixeira et al. 2007). These suspensions 
were heated at 100°C for 30 min in a water bath before total 
RNA extraction was performed. Total RNA was extracted 
using Trizol® standard protocol, based on guanidinium iso- 
thiocyanate and acid—phenol with some modifications (Teix- 
eira et al. 2007). We mixed 200 L of clinical suspension with 
500 wL of Trizol reagent and incubated for 10min at room 
temperature. After addition of 200 wL of chloroform, the re- 
action was mixed vigorously for 10s and centrifuged at 13,000 ¢ 
for 10 min. The upper aqueous phase was mixed with an equal 
volume of cold isopropanol and incubated on ice for 10 min. 
The total RNA precipitate was pelleted by centrifugation at 
13,000 ¢ for 10 min and washed with ethanol. The RNA was 
dissolved in 30uL of diethylpyrocarbonate-treated, sterile, 
double-distilled water and stored at —20°C. 


Reverse transcription—polymerase chain reaction 


The primer combination untranslated region (UTR) 11 —/ 
UTR41+, which corresponds to the conserved region of the 
IBV 3’ UTR and produces a very sensitive reverse transcription— 
polymerase chain reaction (RT-PCR) result (Adzar et al. 1997, 
Culver et al. 2006), was used. Before RT-PCR was performed, 


RNA was heated for 3min at 100°C, followed by 2 min at 
72°C, and finally for 2 min at 50°C. Optimization of a one-step 
RT-PCR requires compromises in the concentrations of PCR 
reagents, particularly the primers used. The RT-PCR was per- 
formed following the instructions in the One-Step Superscript 
II® RT-PCR commercial kit. Briefly, 2 uL of the first-strand 
product was used as a template for amplification in a 50 uL 
reaction mixture containing 20 pmol of primers and 1.6mM 
MgCl. All PCR conditions were the same as described above, 
but the amount of RT-Tag DNA polymerase platinum enzymes 
was calculated to 2.5 units/50 wL reaction mixture. A total of 
10 uL of PCR products was electrophoresed at 100V for 1h 
using a 1.5% agarose gel in Tris-borate EDTA buffer and vi- 
sualized with ethidium bromide staining and an ultraviolet 
transluminator. Gel images were captured using a Kodak 
DC290 digital camera and ADOBE 6.0 software. 


DNA sequencing and phylogenetic analysis 


The PCR products were purified using a commercial pu- 
tification kit (QlAquick PCR purification kit; Qiagen) ac- 
cording to the manufacturer’s instructions and sequenced 
directly with an ABI 3730xl DNA analyzer (Applied Biosys- 
tems). Nucleotide sequences derived from this study have 
been deposited in the GenBank sequence database under 
the following accession numbers: GU323504, Agamia/BR/ 
UNESP/2009; GU323505, APatagonius/BR/UNESP/2009; 
GU323506, Pavus/BR/UNESP/2009; GU323507, Pavus/2/ 
BR/UNESP/2009; GU323508, Penelope/BR/UNESP/2009; 


CORONAVIRUS-LIKE RNA AMONG BIRDS IN A ZOO 


GU323504 (Agamia/BR/UNESP/2009)" 
Gu323510 (Martinica/BR/UNESP/2009)* 


= GuU323509 (Martinica/BR/UNESP/2009)* 
GuU323511 (Rupornis/BR/UNESP/2009)* 
64 GU323508 (Penelope/BR/UNESP/2009)* 
AJ278335 (IBV D207) 


GU323506 (Pavus/BR/UNESP/2009)* 
60) GU323505 (APatagonius/BR/UNESP/2009)’ 


| AJ278338 (IBV HVI-140) 
| AJ278337 (IBV HV10) 


| AJ310642 (TCoV turkey/UK/41 2/00)* 
AJ619580 (Ph/UK/602/95) 


ES AJ619579 _ (Ph/UK/438/94) 
GU323507 (Pavus/2/BR/UNESP/2009)* 


———— 
0.01 
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FIG. 1. Minimum-evolution tree of coronaviruses (CoVs) based on a 146-nt fragment of the 3’ untranslated region from 
eroup 3 CoVs. Evolutionary distances were computed using the Tamura-Nei method and are expressed as the number of 
base substitutions per site. CoVs detected in captive wild birds by this study are denoted with an asterisk. Previously 
published CoV sequences from different sources were included for comparative purposes. GenBank accession numbers are 
shown in brackets. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 
replicates) is shown next to the branches. Phylogenetic analyses were conducted in MEGA4. The scale bar indicates nucle- 


otide substitutions per site. 


ofeese] eesa)s cca] eves daceal sesvadaced | esssdaase al 


5 15 25 35 45 55 
ph /UK/ 602/95 TAATCTGTCT ACTTAGTAGC CTGGAAACGA ACGGTAGACC CTIAGATTTT AATTTAGTTT 
ph/UK/438/94 .... = nauee 2 Slee: mleceanalula im ideeirneaRans mum enatalnalwen dee mmantsvs 
TC oV/UK/ 412/00----.... oo VCE TIS REN WES GRRE ENENN. HSRURE Ee Seeuteweies 
IBV HVI-140 ...... nen ey apaamiansanse fa celeomtnttine wiarsald eamadaiepaldnd - as saieuneiiestane Wasatiabeamarats 
a ee ct SRRRNMNKON REE Re EON ~ aenaeNtehen ara mlNine wares 
IBV-D207 ee rae CSTE eae BSCS SE Keep eeS Beasts: eeaeieeex 
Agamia/BR/UHESP/09 .... jj = «eeevues (3 CAS RMS aaa daa ate gi ears 
APatagmiafiesURese 0S wae SNES ASS See Pees bee 
Pavos/BR/UNESP/09...... Hwee AS SRA saaaans 
Pavos/2/BR/UNESP/09 ... Me kes Tere err eae: 6 Serer ee veCses Kee wee 
PRNRSORE FOOSE 888 cue cai? GaamdeieAirinn «Riel ORR MEME Kiera Re eOee \kiedaeonerarenes 
Martinica/BR/UNESP/09 | = ~~ ........ ce vuhraaees sires caeMenEens | ataniteialee) eaten 
Martinica/BR/UNESP/09 . Ech oes cate we MD Sk Ree Mele ie, Bee es ae 
Rupornis /BR/UNE SP /09 RCH ERR RW: Rea - aeataeeaaae 

beneath aww oul ‘Rews| svediwen') «cesniewedl! neenhecss|| 

65 75 85 95 105 115 
ph /UK/ 602/95 AATTT--GGT TTIAGTTTAAG TTAGTTT-AG AGTAGGTTTA AAGATGCCAG TGCCGAGGCC 
ph /UK/438/94 ..... PTAs: Svineesvas weap ees is eae. Diver gectuawede peeeiuaey , 
TC oV/UK/ 412/00 PR sc ecahniemris “armas Pee acute Din <ctechcaemnie: ax imc cal-ese 
IBV HVI-140 ..... Ts! scat wwe ee Ais soncieaawa) sens G.. 
IBV HVI10 jaws TEA: S4enaeines Gatecun 34 biases skennenawa. wea G.. 
IBV-D207 ~—... Wiha waddadaw Mantas aia interac Bis candles et Seer 
Ag amia/BR/UNESP /09 TUS 4s ne wees Mecéss ie eke was A.. HAS. Kaa G.. 
AP atagonia/BR/UNESP/O9.. .TTA.. cue cee eee wee ee eee ewe ee he SS ee ee ee G... 
Pavos/BR/UNESP/ 09 MTR Beets Seedy wis pet ee Bis poe Megat ene G.. 
Pavos/2/BR/UNESP/09 .G.. **A.. cece Be cee eee Se | Oe Ges 
Penelope /BR/UHE SP /09 eee ccuws eG USGA Bo: civdivases sees G.. 
Martinica/BR/UNESP /09 Ta staaendads areata power aens A.. ere G.. 
Martinica/BR/UNESP /09 PTR ected MO eran sata aNe A.. Wisin <i G.. 
Rupornis /BR/UNE SP /09 TTA Seep ees Soci wees A.. Siivorse® Gains G.. 


FIG. 2. A comparison of part of the 3’ untranslated region from three CoVs isolated from chickens (prefixed IBV), one CoV 
from turkeys (prefixed TCoV), and two CoVs from pheasants (prefixed ph), with the same region identified in the eight CoVs 
isolated from the captive birds described in this study. Dots indicate that all the viruses had the same nucleotide at that 
position. Deletions are marked with an asterisk. Yellow shading indicates the mutations observed in the Pavus/2/BR/ 


UNESP/09 sequence at a given position. 
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GU323510, Martinica/BR/UNESP/2009; =GU323509, 
Martinica/BR/UNESP/2009. The CoV sequences AJ619580 
(pheasant CoV ph/UK/602/95), AJ619579 (pheasant CoV 
ph/UK/438/94), AJ310642 (Turkey CoV turkey/UK/412/ 
00), AJ278338 (IBV HVI-140), AJ278337 (IBV HV10), and 
AJ278335 (IBV D207) were obtained from GenBank for use in 
the phylogenetic analyses. Sequences were assembled using 
SeqMan Pro., and genome annotation was conducted with 
SeqBuilder (DNASTAR, Inc., v.8.0.2). Nucleotide and de- 
duced amino acid sequences were aligned using Clustal-W in 
the MegAlign program (DNASTAR, Inc.). A phylogenetic tree 
was constructed based on amino acid alignments using 
BLASTN analysis (www.blast.ncib.nlm.nih.gov /Blast.cgi). 


Results and Discussion 


Captive animals in zoos are part of a complex ecosystem 
that includes free-roaming animals and their associated in- 
fectious diseases, including respiratory viruses. These results, 
obtained from captive birds in Brazil, represent the first report 
of a circulating CoV inside a zoological park. The CoV RNA 
was detected in eight individual fecal samples (Table 1), 
which corresponds to an individual prevalence estimate of 
1.12% (95% confidence interval 0.5-2.7). Of the samples 
considered positive, two were from black shouldered pea- 
fowl (Pavus cristatus), two were from American purple galli- 
nules (Porhyrula martinica), and one each came from a 
Chestnut-bellied heron (Agamia agamia), an Emperor penguin 
(Aptenodytes patagonicus), a Rusted margined guan (Penelope 
supercilliaris), and a Roadside hawk (Rupornis magnirostris). 
Interestingly, all of the positive samples were obtained from 
adult birds. Phylogenetic analyses were based on a sequence 
of 146nt after the removal of primer sites (Fig. 1). Most re- 
markably, nine mutations were revealed in the Pavus/2/BR/ 
2009 (GU323507) sequence, despite many reports that con- 
sider the 3’ UTR region to be highly conserved among group 3 
CoVs (Cavanagh et al. 2002). One hypothesis that could ex- 
plain these mutations is the high density of different species of 
birds located in the same geographical area, which allows the 
virus to circulate among species (Jackwood et al. 2010). On the 
other hand, the natural habitat provides more distance be- 
tween species, particularly between predators and prey. In 
addition, the sequences of the CoVs isolated from the two 
peafowl samples clustered separately from the British de- 
posited sequences, as illustrated in Figure 1 (Cavanagh et al. 
2002). Sequences within this cluster were relatively homoge- 
nous between two European samples (Fig. 1), in accordance 
with previous studies (Cavanagh et al. 2002). 

Samples obtained from the Chestnut-bellied heron, two 
American purple gallinules, a Roadside hawk, and a rusted 
margined guan clustered together in a single clade, close to 
the sequence of the IBV D 207 strain (Figs. 1 and 2). In the 
present study of samples collected from 25 different species 
(kept in isolated locations within the zoological park), CoV 
RNA was detected in two Galliformes, two Ciconiiformes, 
one Gruiforme, and one Sphenisciforme. This demonstrates 
the high diversity in the susceptible species, as described 
previously (Woo et al. 2008). In fact, CoVs have been also 
detected in wildfowl species, rock doves, wild peafowl, 
waders, Amazon parrots, and some passerine species (Cava- 
nagh 2005, Liu et al. 2005, Gough et al. 2006, Woo et al. 2008, 
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Hughes et al. 2009). The hypothesis that wild free-flying birds 
present in zoos may often act as vectors flying between en- 
closures that house the exotic birds and transmitting infec- 
tions is not fully understood. 

To understand the potential role of captive populations as 
reservoirs for CoVs, further studies are needed not only in 
South American zoo parks, but also worldwide. This is the 
first description of the molecular detection of group 3 CoVs 
related to captive birds in South America. The genetic diver- 
sity found in the black shouldered peafowl samples raises 
interesting questions about the 3’ UTR conserved region, its 
evolution, and viral circulation among different species inside 
ZOO parks. 
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